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RANDOM EVEN IB WERE. SIMULATED BY COMPUTER AND SUBJECTED TO VARIOUS 
STATISTICAL METHODS TO EXTRACT IMPORTANT PARAMETERS. VARIOUS 
I OUTS OF CURVE FITTING WERE EXPLORED, SUCH AS LEAST SQUARES, 

LEAST DISTANCE FROM A LINE, MAXIMUM LIT FLIHOOD. PROBLEMS 
CONSIDERED WERE DEAD TIME, EXPONENTIAL DECAY, AND SPECTRUM 
EXTRACTION FROM COSMIC RAY DATA USING BINNED DATA AMD DATA FROM 
INDIVIDUAL EVENTS. COMPUTER PROGRAMS, MOSTLY OF AN ITERATIVE 
NATURE, WERE DEVELOPED TO DO THESE SIMULATIONS AND EXTRACTIONS 
AMD ARE PARI I ALLY LISTED AS APPENDICES. THE MATHEMATICAL BASIS 
F-OR THE COMPUTER PROGRAMS IS GIVEN IN THE TEXT OF I HE REPORT. 
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2. RELATIONSHIP OF BINOMIAL, POISSON, AND NORMAL DISTRIBUTIONS 

IF WE DIVIDE A TIME INTERVAL T INTO n EQUAL SMALL 
SUB INTERVALS OF LENGTH Dt , A RANDOM EVENT IS EQUALLY LIT ELY TO 
OCCUR DURING ANY SUB INTERVAL OF DURATION Dt. THE PROBABILITY Dp 
OF AN EVENT DURING AN ( GIVEN SUB INTERVAL Dt IS GIVEN BY THE 
LX PRESS I ON, 


Dp - L Dt 


2 . 1 


WHERE L [3 THE 
CO SMALL THAI 
UNLII ELY. THE 
GIVEN iiMILRVAL 


AVERAGE NUMBER OF EVENTS PER UNIT 
TWO COUNTS DUPING TIC- INTERVAL Dt 
PROBABILITY THAT AN EVENT WILL NOT 
D i_ IS GIVI.N BY I HE EXPRESSION, 


TIME AND 
IS HIGHLY 
OCCUR HI 


DT IS 


ANY 


Dp = i 


DT 


rill. F RUBAI ITY I ‘ ( I i , I ) 01 
DfCfRj mn I UN, 


N EVEN IS IN f I f !F 1 


EG t) EVEN BY IT IF UFIIUMIAI 


PEN, I) - I l ft-DL T N * T ( J - E) J . I (n N>T« n'/T( n-N) 1 N'l - 

1 1 1 FQ AND 101.1 OWING, 1 1 Ilf SYMBOL * DEMO E f:S Mill. T ET'L I CAT TON AMD 

1 1 IF SYMBOL I ’R! Cl Dl S AM i XPONI iM I . P(N,'i> IS IMF Nl h TFKH OF TUP 
PCI VN0MIA1 IMfil i;j I ft r>1 ) I n THIS DfSlRlBUilON IS Ei-IERLI ORE 
NORMA! I /if), AND AN/ R EG / K I RU I I ON D! LIVED T ROM IT TN A E Til FT E NO 

process siioui d be nof<itai f/fd. 

1111 POIGGOil DESILTHH10N IS ! "AS 1 1 Y f)B i A f PIL'D T F^OM M IE 
BI NOT! I i'il DfSERlBUl I ON BY USEMO ! HE S f ERL tf 10 APPROXIMATION, 

100 n 1 -- (1/2) IDS <2 P \> i- < n i 1/2 ) LOG ri - n , 2.4 

AND AIT POX 1 MAI IONS I HE 

100 (n-N) - 100 li.(1 I! / n ) I - I 00 n - n/N . 2.5 


TliESi: APPROX l MAT EONS ARE QUITE OOCJD WHEN n , N. 
IMA I URAL I HOAR I HIM Of EU 2. .1 AND USE. EQUATIONS 2.4 
OB I A 3 T! El IE PC! 1 3 3 ON D I S I R I BUI 1 ON , 


IF WE TAI E THE:' 
AND 2.5, WE 


1 1 ( M , T ) •= I. ( L T ) NTn-L EIX'P ( Li-T) J/N 1 


6 


HIE AVERAGE NUMBER OF EVENTS N IN TIME T IS ( L s T ) AND THE 
STANDARD DEVIATION FOR THE POISSON DISTRIBUTION IS SQR ( L*T ). 

AS THE NUMBER OF EVENTS N BECOMES L.ARGE , THE POISSON 
BECOMES EQUIVALENT TO THE NORMAL DISTRIBUTION. IF WE TAKE THE 
NATURAL. LOGARITHM OF THE POISSON DISTRIBUTION AND USE THE 
STIRLING APPROXIMATION AND THE APPROXIMATION, 

LOG C 1+X/ (L«T) 3 = X / L T - C(X/L*T)'23 / 2 , 2.7 


WHERE 
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X 


r-4 - l i 


1 


2.8 


WE FIND THAT 

P ( X , T > ■- SOP C 1 / (2 PI Lit T 1 k EXP C- X 2/ ( 2 L*7 ) 3, 1.9 

WHICH IS FHE NORMAL D [ ST R I BUT i(3N WITH AVERAGE LT ( X = O ) AND 
STANDARD DEVIATION 30R ( L T > . T A RLE 2.1 COMPARES THE DINOMIAL , 

PO J BOON AMD NORMAL DISTRIBUTIONS FOR N ■= 30 , T -= 5 , AND L - 3 . 

3 Dl AD TIME CURE!" CT IONS FOR RANDOM EVENTS. 

SUPPOSE HIAf AN EVENT DETECTOR (COUNTER) HAS A DEAD TIME t 
AF ! ER I ALT I COUNT. IF M ' RANDOM LVLNT3 ARE OBSERVED DURING TIME T, 

WE i NOW I HA f THE GOWN I ER WAS DEAD ( J NOFT.RA T T VE> FUR A f THE M ' H . 

1 HI AVERAGE Ni JITBER OF" EVENTS EXPECTED DUR [NO THIS DEAD TIME WOULD 
BE LkW'H. AMD WOULD Ai SO OBEY A POISSON DISTRIBUTION. THE NUMBER 
OE E VEIT IS M FS E Si IMA I FID 10 BF GTVETI BY THE EQUATION, 

|\j - j,| - i I. N ' L . 3.1 

IT DUI S NO I MATTER MOW OF I E'N OUR COLJNFER IS FUF<NI£D ON OR OFF. IF 
THE EVENTS ARE RANDOM AND L IS A C0N3TAN F , ED 3.1 IS VALID WHEN 
N REPRESENTS AIL I HE OBSERVED COUNTS AMD I REPRESENTS FHE TOTAL 
f THE TilAf THE COUNTER JG ON. J E WE MAI E IF IE ASOUN T I ON THAT N 3S 
F 01 LAI 10 List „ LU 3.1 BE COM! 'S 

LJ«T - N't L * N ' k | 3.2 


OR 


L-N'/tTN'tt) . 3.3 

EO 3.3 IS THE SIMPLE STATEMENT THAT HIE BEST ESriMATE" OF HIE 
MATURAI COUNT RATI L IS 1 1 IE IOTA!. NUMBER Of COUNTS DIVIDED BY THE 
TOTAL TIME THE COUNTER IS ON. IT IS POOR PRACTICE TO REPLACE L IN 
EO 3.1 BY T! / T AND USE THE RESULT, 

N ~ N ' / ( 1 - N't/ T ) , 3.4 

AS A CORRECTED COUNT. THIS REACT ICE OVERCORRECTS WHEN N ' IS 
LARGE AND UNDERCORRECTS WHEN N ' IS SMALL. 

THERE ARE A NUMBER OF TECHNTOUES WHICH MAY BE USED TO STUDY 
COUNTER DEAD TIME. A FIRST APPROACH IS TO OVERWHELM THE COUNTER 
WITH MANY MORE EVENTS THAN IT CAN POSSIBLY COUNT. IF THE EVENT 
RATE L IS RELATED FO THE COUNT RATE R BY TFTE EQUATION, 

P ■= L / ( 1 4 L*t ) , 3.5 

WHERE t IS THE DEAD TIME AFTER EACH EVENT, THE MAXIMUM POSSIBLE 
COUNT RATE R (MAX ) IS SEEN TO BE GIVEN AS 1 / t , SO THAT THE DEAD 
TIME t = 1 / R (MAX ) . IF THE COUNTER DETECTS RADIATION, A MORE 
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TABLE £.1 


A comparison of binomial, Poisson, and Gaussian 
distributions for n = 30, T = 5, and X = 3 
corresponding to an average count of 15 


N 

B(N) 

P(N) 

G(N) 

1.000 

0.000 

0.000 

0.000 

2.000 

0.000 

0.000 

0.000 

3.000 

0.000 

0.000 

0.001 

4 .000 

0.000 

0.001 

0.002 

5.000 

0.000 

0.002 

0.004 

6.000 

0.001 

0.005 

0.007 

7.000 

0.002 

0.010 

0.012 

8 .000 

0.005 

0.019 

0.020 

9.000 

0.013 

0.032 

0.031 

10.000 

0.028 

0.049 

0.045 

11.000 

0.051 

0.066 

0.060 

12.000 

0.081 

0.083 

0.076 

13.000 

0.112 

0.096 

0.090 

14 .000 

0.135 

0.102 

0.100 

15.000 

0.144 

0.102 

0.103 

16 .000 

0.135 

0.096 

0.100 

17 .000 

0.112 

0.085 

0.090 

18.000 

0.081 

0.071 

0.076 

19.000 

0.051 

0.056 

0.060 

20.000 

0.028 

0.042 

0.045 

21.000 

0.013 

0.030 

0.031 

22.000 

0.005 

0.020 

0.020 

23.000 

0.002 

0.013 

0.012 

24 .000 

0.001 

0.008 

0.007 

25.000 

0.000 

0.005 

0.004 

26.000 

0.000 

0.003 

0.002 

27 .000 

0.000 

0.002 

0.001 

28 .000 

0.000 

0.001 

0.000 

29.000 

0.000 

0.000 

0.000 

30.000 

0.000 

0.000 

0.000 
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SUITABLE APPROACH INITIALLY IS TO USE A RADIATION SOURCE AND THE 
INVERSE GOUARC LAW LOR CAL ] DRAT ION „ THE EVENT RATE AT THE COUNI ER 
JS GIVEN AS 


WHERE T tS PROPORTIONAL TO THE SOURCE INTENSITY. FROM EOS G. 5 AND 

« < » 1 


R =• L / ( 1 " LkT ) - 1 / [ ( r 2 / E > i- t J , 


2. 7 


GO I HAT 


- - R*l *t 


9 


A PLOT OF R*r 0 VERSUS R SHOULD BE A STRAIGHT LINE WITH SLOPE 
I * AND fNTERCLPIS, 


R *r 2 - I ; R --= I / t . 'LJo 

IE THIS IS NOT SO, DEAD TINE t IS NOT INDEPENDENT OF COUNT RAIL. 

4. COMMON ST A I 1ST J CS FOR EXPONENTIAL DECAY 

THE PROBABILITY Dp OF A RAD TO ISOTOPE DECAY DURING THE TIME 
INTERVAL D i IS GIVEN BY THE EUUAIION, 

Dp I.HNkDT , 4.. I 

WHIRL N TS HIE NUMBER OF A I DNS AVAIL ABIE FOR DECAY AND L. IS THE 
DECAY CONSTANT. IF WE DIVIDE TIME I INTO n EDUAl. INTERVALS DT, 
r I IF" PROBABILITY OF NO DECAY IN TIME T FOLLOWED BY A DECAY IN TIME 

nr ro given as 


P <r>*I>t -II- I. k N fc-DT ) n kL*N*DT . 4.2 

AS n BECOMES LARGE THIS REDUCES TO 

PIT) fc-riT = L EXP ( L x-TTr- I > TM_>N*dT . 4.G 

THIS DISTRIBUTION OF THE TIME BETWEEN DECAYS HAS AN AVERAGE OF 

I *N AND A STANDARD DEVIATION OF L»N. GIVEN THE TIMES T(I) 

ASSOCIATED WI TH A NUMBER OF DECAYS N(I), A NUMBER OF TECHNIQUES 
HAVE BEEN USED ID COMPUTE THE ORIGINAL NUMBER OF ATOMS N 1 0 ) AND THE 
DECAY CONSTANT L. IF THERE ARE MANY DECAYS PER UNIT TIME, IT IS 
CONVENIENT TO WRITE THE DIFFERENTIAL EQUATION, 


dN / cl T •= - L N , 

WHICH HAS THE GOLUFION, 

N(T) = N<0> EXP ( - L T ) 


4.4 


4.5 
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SINCE THE COUNT RATE IS PROPORTIONAL TO N (T) , A PLOT OF LOG 
COUNT RATE VERSUS TIME HAS A SLOPE OF -L AND AN INTERCEPT L*N(0>. 

IF EVENTS ARE FEW, WE NEED TO DEVELOP TECHNIQUES TO EXTRACT 
AS MUCH INFORMATION AS POSSIBLE FROM THE DATA AT HAND. TO TEST 
THE VARIOUS TECHNIQUES, DECAY DATA WAS GENERATED BY COMPUTER AND 
THE DATA WAS PROCESSED BY EACH TECHNIQUE TO SEE HOW WELL IT WOULD 
DO. 


A COMPUTER GENERATES RANDOM NUMBERS X ( I ) DISTRIBUTED 
UNIFORMLY IN THE INTERVAL BETWEEN ZERO AND UNITY. TO PRODUCE 
NUMBERS CORRESPONDING TO RANDOM DECAY TIMES, WE EQUATE THE 
DISTRIBUTIONS, 


P(T) dT - F(X) dX, 4.6 

WHERE P'D COMES FROM EO 4.3 AND F(X> IS UNITY. THE INDEFINITE 
INTEGRAL 01- BOTH SIDES OF EO 4.6 YIELDS 

T ~ - ( 1 / N « L ) LOG ( 1 - X ) . 4.7 

BECAUSE X IS A RANDOM VARIABLE, WE CAN WRITE 

1(1) - -- ( 1 / N < 1 > *L ) LOG ( 1 - X(I> ) , 4.0 

WHERE 


NCI) - N ( t"! ) - T +• ! _ 4 _ 9 

JF EVEN IS ARE FEW, WE CAN RECORD THE TIME T(I) FOR EACH DECAY AND 
FT I i .0 4.5 TO THE DATA. THE RESULTS ARE GENERALLY POOR IF WE USE 
LEAS I SQUARES FITS. USING THE DISTRIBUTION OF ED -4. 3, WE CAN 
GAL CULAIE THE AVERAGE f I ML t(I> BE TWEEN DECAYS TO BE 

t ( i ) = L / r. ( N ( 0 > - ] + 1)S*L I . 4.10 

AT I ER REARRANGING, 

IM-F1) ---- ( N ( O ) + 1 > (!) -- 1 / L . 4.11 

A LEAST SQUARES CURVE FIT OF 1 p t ( I > VERSUS t(I> SHOULD YIELD A 
STORE UP M (U)il AND AN INTERCEPT OF -1 / L . THIS OFFERS LITTLE 
IMPROVEMENT. WE COULD REARRANGE ED 4.5 TO THE FORM, 

T<3> ~ (J/E) LOO N(D> - ( l / L > LOG N ( I > . 4.12 

A LEAST SQUARES FI I UF TCI) VERSUS LOG N(I> YIELDS A SLOPE OF 
-i/L AND AN INTERCEPT OF (1/L) LOG N(O). 

WE MIGHT THINK THAT THE METHOD OF LEAST SQUARES GIVES TOO 
MUCH WEIGHT TU ERR AT I C (RARE) EVENTS AND THAT SOME SCHEME TO GIVE 
LESS WEIGHT 10 POINTS DISTANT FROM THE CURVE FIT MIGHT BE BETTER. 
A SQUARE ROOT AND FOURTH ROOT OF THE SUM OF SQUARES RESULTED IN 
SMAL L T MT ‘ROVE ME! NT S . 
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UGi.D TO DERIVE I- AND N 
i Vf m is, A>i SCR I Hi" D r:Y I 
] i , NAG A l--RfJL.r.r« J I. ITY 

such a ■>(• on: nci can hi 


Hr, THE ME 1 NOD 
(!>•. n If, AEG 
A! l!) i ( l > , I Hr 
CLOG!" IN lilL- 


nr ill' ,x ] mi iii i.nri moon wag 

LINED I NAT A Sl-'OUENLE Ur DK CAY 
! j I ill £! L I W; I N LVQN G 1 AND 
!'IAX 1 MNi'l. NIL PROBABILITY FOR 


WR I I If M AG 


WHI 


i 1 1 
1 


rr t n ) , t f ?> , „ . 'j - ( i > vp (?) - r ( 3> * . . 

f I'G A CONST ANT Oi I ‘F-TJI OR P I (JI4ALTTV 
P (I ) --- I ; LNM ( i ) ft EXP [ N(l)M ittG) 3 


AT ID 


4.. 1 3 


4 . J 4 


'NIL'-. l-*r<f );>( J.' ; I 3 Mf JIM- NAMAGAHI t 


A I GOAR I 1 f IN 


lug i- n ( i ) , i- < > 


3 — Ml IN ( f) N DHL a *1 s» N ( 1 > 


iJGNin AN I 11 LA I I AN GLUE- 1 II- W! 

max rnizr iiiN ruNcrjoM. 


Di IOOGP 


GI 


IN INF FORM, 

I N < f ) M. » I < f > < 4. I’G 
I' NO) AND I WIIIGII 


LLGUI I ' > 1 )G I l\K / 
IT IG SEEN FNA f f 1 113 
GUI T R3 OR I OR N (O) 

rs nor i i iwn il e lf tni 

IAI'.iI GNOWG Nlh 


1 1 if. vai [ ( iur> m:i in j lm h r. ari , ;i n iwn 
in : n n id or maximum lii i i mono ig 


Cl I 1 AH! L 4 - l . 

cgng cut run y 


!“« 1 'I AC NIG I LOG I . II ‘HAG I G ON DEGfANT NOENIG 
1 1! rilOD GI' NAXCNuH i It- r I r 1100.0 eg ngi ljgfd,, 
RrGOLNG ill NGINC INI 1> 1 1 i> Ed N 1 ML1 HODS ON NIC 


GAI il. 


GI DA I A. 


mG f.N AG ! BE , II IG NON h 1 1 In I 10 4.10 CONI D 
PG 3 0001 1 D 1 0 I P I ON I J ON IP I. * N C T ) f S I v'i. PLACED HY I .. 
f. LVI 1 1 AG 


Hh OGf D I OR A 
f! IF RFC GUI I IG 


I 00 IT t ( I > , t (?) , . . I = GUM (1 ) I 00 I'M i . ft l ( I ) 

- No LOG I x|. -I M 


4 


1 h 


WlirKI I M> II IE 
riJNCT LON OCCORG 
BE I WEEN L VENTS 
DIMLY INC lOTAl. 


101 AL Til-ir FOEv N COUNT G. 1 1 111 MAX EH I7A I LON OF" THEG 
WNL'N L-N/T, SO TIIAE THE: DIG FR I DOT rON 0I : TINE 
I G ME ) 1 I KOI" I "IJL IN G I UDY I NO A POM ISDN D E G I H 1 HU 1 I ON . 
FI Ml I AND THE: IOIAL COON Ff» N 3S IMPORTANT. 


TEC A 1 ML N I 0! " PE JWER GI T- L 1 RUM DA I A 


EVEN EG DULY ENG TNI". DIG FR I HUT EON, 

P ffl) JvtiK = A x F-I iM HE. , 5.1 

ARE Cl -I ARAL TE.R 1ST IC OF COSMIC RAYS AMD RAD TATI ON BELTS, WHERE E I 
TOE PARTICLE ENERGY AND P(E>#<iE IS THE PROBABILITY OF OBSERVING 
F'ARTICLE WITH ENERGY BETWEEN E AND Ei-dE. f HE NORMALIZATION OF 
T FIT S FUNCTION RESULTS IN I TIE LOUALITY, 

A - ( B - 1 ) / ( LL ( BH) -UL ( B+l)) , 5.2 

WHERE LL IS THE LOWER ENERGY LIMIT AND UL IS THE UPPER ENERGY 
LIMIT. USING THE METHOD OF SECTION 4, TINS SPECTRUM IS S EMULATED 
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TABLE 4. 1 


ORIGINAL PAGE IS 
OF poor Quality 


SUMMARY OF FIND I MGS ON CURVE FITTING 


N<0> =• 

5 

LAMBDA 

= 1 

SEQUENCES =• 7 


METHOD* 


LAMBDA 

S. D. 

M (0) 

S. D 

i 


2.45 

cr *-y~r 

3.68 

1 . 14 

n 


0.71 

0. 5S 

5. 46 

1 . 62 



l . 0 1 

0 . 62 

6 . 13 

2.66 

4 


1 . 24 

3 . OS 

7. 49 

6.61 

N ( t > > =• 

1 0 

LAMBDA 

- 1 

SEQUENCES •= 7 


METHOD* 


L AMBDA 

S. D. 

N ( 0 ) 

S. D 

1 


2 . 03 

3 . 03 

0.56 

0. 69 

n 


0. 93 

0 . 50 

3 0. 65 

'“i '“I -*7 

“•* 


O. 97 

0 . 52 

13.32 

2 . 50 

4 


0.97 

u. 18 

10.46 

1 . 19 

5 


l.i>4 

0 . 28 

1 > ). 06 

1 . 96 

N (O) - 

30 

LAMBDA 

= l 

SEQUENCES - 7 


mi: i i-idd* 


LAMBDA 

G. D. 

N ( 0 ) 

S. D 

i 


2 . 0 

0.45 

27 . 92 

1 . 44 

r-i 


' 73 

O . O 7 

nr) ^ 

4 . 05 

Z 


0.99 

O.'-'O 

29 . 65 

4 . 30 

4 


1 . ! >0 

0 . 1 ? 

30. 1 3 

0 . 29 

N(0) - 

So 

LAMBDA 

-- 3 

SEQUENCES 7 


METHOD* 


LAMBDA 

n. n. 

H(o) 

S. D 

3 


t . 39 

0 . 22 

47.02 

3 . 67 

' > 


0 . 9 1 

O.J4 

4 7 . 1 2 

3. O L 



0. 93 

0.14 

40.06 

2.87 

4 


i . Of > 

0. 1 3 

5c. 51 

0.03 

15 


L . 0 1 

0.33 

50 .46 

1 „ t >0 

6 






METHODS 





' 

1. L.S.F. 

ru 

IM'AII(J) - 

( N ( 0 > +• 

i MTAU(I) - i/L 


2 . L.S.F. 

T 0 

LOG ( N((>) 

- I ) — - 

L * T ( I ) 4 LOG ( N ( 0 ) ) 


3. L.S.F. 

10 

1(1) -- (-1/L)*L0G ( 

N(0> - 1 ) h (3/L) LOG 

(N (U) 


4. MAXIMUM L. II EL I HOOD ITERATION 

5 . LEAST D ] S 1 ANCE 1 T ERA F 1 ON 

6 . ROOT MEAN SQUARE ITERATION 

L.S.F. - LEAST SQUARES FIT ; L - DECAY RATE CONSTANT 
TAU(I) ~ TIME BETWEEN EVENTS 1-1 AND I ; N(0> - ORIGINAL 
NUMBER OF EVENTS ; T < I ) - TIME WHEN EVENT I OCCURS. 


XXXVI I- 7 



TABLE 4 


OiT .. 
2 'K' 


COMPARISON Of SOME METHODS USING THE SAME SET OF DATA 


SEED 

= 

11 1 1 1 EVENTS ■= 30 

LAMBDA = 1 

N (0) 


LAMBDA 

METHOD h 

EDO AT iOLHt* 

OS. 9 


1 . 27 

1 

1 

27.5 


O. 88 

1 

2 

28. 8 


0. 93 

1 

r t 

29. 3 


2 . Ob 

n 

1 

29 . 7 


0. 9 7 

r t 


29.8 


0. 90 

n 

“•j 

29. 1 


2. 99 


1 

23.9 


O. 96 

-*• 

**b 

3 3.9 


1.17 

“1 

j 

29. 2 


1 . 99 

4 

3 

2 7 . b 


0. 88 

4 

* b 

30. 6 


1 . t >1 

4 


29 . G 


1.14 

!j 

4 

ME i HODS 





I.E AST SUM 

or 

COUARFS 

OF (Y( D -N k X ( D LO 

LEASE SUM 

OF 

ABO (YU 

) M#X ( I ! B) 


LEASI SUM 

or 

SUR I ABS 

( Y (1 ) -M* X (1 ) 

+ B ) ) 

LEAST SUM 

or 

DISTANCES OF POINTS 

FROM LINE 

MAXIMUM 1. 

U L 

L I HOOD 



r.OUAl 10N3 





1 1 1 ( J ) 

(M (o) + J > *t (I) -t/L 



2. LOO N ( I ) =• ~L*T(J> + LOG N(0> 

3. T (1 > - -I. LOG MCI) J/L -i r LOG N(o> I/L 

N ( J ) « t VENTS REMAINING 

1(1) == rillf WHEN Hit EVENT OCCURS 

L(l) •--- T ( 1 ) -T ( E - J ) 
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BY THE RANDOM NUMBER STATEMENT, 

E ( I ) ■= < LL' (— B+l ) - C (Bd)*Xd)/AD )'(1/(B-1)) , 5.3 

WHERE X(I) IS A RANDOM NUMBER DISTRIBUTED UNIFORMLY IN THE 
INTERVAL CO, ID. THE MAXIMUM LI I EL I HOOD STATEMENT IS WRITTEN AS 

LOG P ( 1 , 2 , . . > = SUM ( I ) C LOG A - B*LOG E(I) . 5.4 

TABLE 5.1 SHOWS THE EFFECTIVENESS OF THE MAXIMUM LIVELIHOOD 
METHOD OF EXTRACTING THE EXPONENT B FOR THE DISTRIBUTION. 

IF DAI A IS BINNED, THE PROEABILTY P f E ( I > , E ( I + 1 ) D OF AN AN 
OBSERVED PARTICLE HAVING ENERGY IN THE "BIN" BETWEEN E(I) AND 
E d H ) IS NIL ] NT F- ORAL OF THE NORMAL I Z ED D I STR JBLJT I ON FUNCT 1 ON 
BETWEEN THESE ENERGIES OR 

PLEd ) „Ed i 1 ) D = (A/(B 1 ) > CE ( I ) ( - B M ) Ed + 1) < -Eh 1 > D . 5.5 

THE MAXIMUM Ltm.fllDOD SfATEMENT FOR A DINNED DATA IS 

10(5 I-' I J ,2, S, „ . . 1 - SUMd) Bd) LOG i 0 3 5.6 


Wl IF.RE 


0=fC E(I> ( B+D-Edd) ( - B+l ) 1 / C LL (-B-d)-UL (-Bd) 3 . 5.7 

THE RESULTS OF SOME COMPUTER EXTRACTION OF EXPONENTS FOR POWER 
LAW SPEC IRA ARE GIVEN IN TABLE 5.1. IT SHOULD BE NOTED THAT 
BINNING IS A WAY IT INTRODUCING UNCERTAIN TIES IN THE DATA OR A 
WAY OF THROWING AWAY INFORMATION EITHER BECAUSE OF CONVENIENCE OR 
NECESSITY. HIE RESULTS OF TABLE 5.2 ARE FOR BINS OF EQUAL SIZE. 
THE FNERGH-S Ed) DEFINING BINS COULD HAVE BEEN CHOSEN BY 
GEOMETRIC PROGRESSION OR INTERVALS OF EQUAL PROBABILITY OR ANY 
OTHER WAY DESIRED. 

AS A SUMMARY, WE CAN SAY THAT DATA SHOULD NOT BE BINNED 
UNLESS IT IS INESCAPABLE. 

6. RANDOM ERROR AND SPECTRUM EXTRACTION 

FOR VARIOUS REASONS, THE RESPONSE OF THE DETECTING 
l NS T RUMEN I S IS DIFFERENT FROM THE REAL ENERGY SPECTRUM . THE 
LEAST TROUBLESOME ERRORS ARE THOSE WHICH RESULT IN AN ENERGY 
UNCERTAIN ! Y OVER A LIMITED RANGE. BINNING, FOR EXAMPLE, 

INTRODUCES AN UNCERTAINTY EQUAL TO THE BIN WIDTH, AND WE CAN 
SIMULATE QUITE EASILY WHAT BINNING DOES TO THE SPECTRUM 
EXTRACTION. GIVEN A BIN DISTRIBUTION AND THE ABILITY TO SIMULATE 
THE EFFECT OF BINNING ON A GIVEN SPECTRUM, WE CAN ALWAYS FIND A 
SPECTRUM WHICH REPRODUCES, WITHIN' THE VAGARIES OF STATISTICS, THE 
BIN DISTRIBUTION. SOME INSTRUMENTS WILL DETECT PARTICLES OF A 
GIVEN ENERGY Ed) AND INDICATE A DISTRIBUTION OF ENERGIES E(J> 
SUCH THAT 0 I E(J) < Ed). 1 HE PROBLEM Or EXTRACTING THE ORIGINAL 
SPECTRUM THEN BECOMES DIFFICULT OR IMPOSSIBLE. 
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TAFJLF. 0.1 


MAXIMUM L If ELI HOOD TREATMFN f OF POWER DiSTRI DUTTONS 


Oi-F D -= 0 1 0 1 12 

OEr-D = 

12 l 12 1 2 

HI- hi) -- 

2 1212 

LI -1 UL-5 

LL- 1 

UL=10O 

LL- : 3 

l il. = l UOU 

F X F ' -= 2.7 

EXP = - 

-2. 7 

EXP - - 

'~V 

/ 


even rs 

EXP 

EVEN f S 

EXP 

EVENTS 

EXT' 

1 1) 

' > • 

J 0 

-2. 7 

lo 

-3.o 

2u 

2. 9 

20 

-3 . 0 

20 

- 3. 1 

'■0 

-2. ! ? 


-2.9 

70 

-2. 9 

40 

2. 7 

550 

H O 
j. « U 

5u 

2. 7 

Uo 

2„ n 

1 oo 

2. £1 

1 U’ 

• 7 ' 

1 U 1 > 

loo 

• \ * 
». ■ / 

2 f )u 

- 2.9 

2 1 u_i 

a ' u £> 



<1 ( a ) 

' > O 

n / 

0* > 

r * V 





TAD! C 

3. 2 


PL 

HUE IS or 

MAXIMUM 

LIT LL THUOL) 

IFF A 1 ML NT 01 HINMfMG 

ERRORS 

EVEN 1 S 

l.l. 

UL 

BINS 

EXP (GIVEN) 

EXP (FIT 


j 


4 

2 .. 7 

3 „ 1 

Ui ) 

1 

"T« 

0 

-2.7 

- 7 . 2 

1 oo 

J 

3 

lo 

-2.. 7 

- 7. 15 

2ot ) 

J 


20 

-2.7 

-2. S 7 

'f/Oi ) 

1 

5 

20 

2. / 

2.82 

300 

1 

b 

TO 

-2. 7 

-2.81 

500 

t 

1 O 

2o 

.. '~l ~7 

- 2 . 79 

1 00 

t 

20 

TO 

-2.7 

'•7 17 ' 7 
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AS A FIRST EXERCISE WE WILL LOOK AT THE EFFECT OF A 
MEASURING UNCERTAINTY WHICH IS PROPORTIONAL TO ENERGY, NAMELY, 

E (F) = E<I>*( 1 -f FM0.5-RND)) , 6.1 

WHERE Ed) IS Ti-IE INCOMING ENERGY, E (F) IS THE READING FROM THE 
INSTRUMENT, F REf-'RESENTS A RELATIVE SPREAD AND RND IS A RANDOM 
NUMBER UNIFORMLY DISTRIBUTED BETWEEN ZERO AND ONE. FOR EXAMPLE, 

IF < = 0.1 THEN 0.75*C(1> E (F) 1.05*E(I>. TABLE 6.1 SHOWS THE 

t. FFECT OF INTRODUCING A RANDOM ERROR SUCH AS INDICATED IN EO 6.1. 
IT IS GEI.N THAT A RATHER LARGE ERROR OF THIS TYPE CAN BE 
T OLE RATED WITHOUT SERIOUS DEGRADATION OF OUR ABILITY TO EXTRACT A 

gpf: c trljm . 

TO ILLtJSIRAFE HIE SECOND F IND 01- ERROR, WE IMAGINE THAT THE 
INCOMING ENERGY E IS DEGRADED TO A READING E' ACCORDING TO THE 
D 1ST RI DU I HIM, 

PIE, El') - A *■ k E (-B) ) *A ' rfl.XP I - (L--E ' ) / <h*E) , 6.2 

WHERE AKE(-B) REPRISE NTS THE INCOMING SPECTRUM AND THE REMAINING 
FAC I OR Rl.TREGL'NT S THE DEGRADATION OF THE INSTRUMENT. THE PRODUCT 
h r E IG THE AVERAGE DEGRADATION C -E ' . THE RESULT OF USING THIS 
DISTRIBUTION IS SHOWN IN FABLE 6.2. IT TS SEEN THAT THIS TYPE OF 
ERROR MAILS THE SPECTRUM LOOT STEEPER BY TRANSFERRING EVENTS FROM 
1 1 1 TIMER I 0 I OWE I < I: TIERS I EG . 

7. COSMIC RAY DATA T ROM CLRENI OV COUNTERS 

WHEN A HIGH ENERGY NUCLEUS WITH SPEED V AND ATOMIC NUMBER Z 
1 RAVE! ROES A REFRACTIVE MEDIUM WITH REFRACTIVE INDEX N, A LIGHT 
PULSE IS GENERA I ED ACCORDING TO THE FORMULA, 

I. -- I *(Z 2>*( 1 - CC/(V*N>> 2 > , 7.1. 

WHERE C IS THE SPEED 01 LIGHT IN VACUUM AND I IS A CONSTANT. IT 
IS SEEN FROM EO 7.1 THAT LI MAX), THE MAXIMUM LI GIFT PULSE 
AVAILABLE WHEN V/C -= i, IS GIVEN AS 

LIMAX) - MKZ 2>/( 1 1/N 2 ) . 7.2 

THE INDEX OF REFRACTION N OF THE GAS IN THE COUNTER WAS 1.00115 
SO THAT L (MAX) / (T * (Z 2)) ■= 435.53 . 

THE CERENF OV COUNTER ACTS TO CONVERT THE ENERGY SPECTRUM, 

N (E) - A * E''(-B) , 7.3 

TO THE PULSE HEIGHT SPECTRUM P(L> WHERE 


N (E) dE •= PIL) dL 


7. 4 
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TABLE 6. 1 


EFFECT OF SIMPLE ENERGY UNCERTAINTY ON SPECTRUM EXTRACTION ► 
P(E> = A E (~B) B - 2.7 LL = 1 UL =5 



E ( F ) = E ( I ) *■ ( 1 i h# ( 

O. 5 

-- RND ) ) 

EVENT S 

ERROR PARAMETER 

h 

EXTRACTED B 

10 

O. 1 


'~t / 

.►_ « O 

10 

0. 2 


2 .. 09 

in 

0. 3 


2.5 7 

i<> 

o. o 


j 1 .. r Ci / 

1 0 

0. 9 


2.63 

1 0 

3 . 9 


3.84 


TABLE 6.2 


ETFECT OF EXPONENTIAL DEGRAD A I ION ON OPEC i RUM EXTRACTION 
P(E> --= A E < --B ) B = 2.7 LL. = 1 UL ■= 5 
P(E.,F'> - A*E (-R)-A’EXP C- <E- E ' > / <h*E > 3 


EVEN TO 


ERROR PARAMETER h EXTRACTED B 


LO 

O 

10 

0 

10 

0 

1" 

0 

1 0 

0 

7i> 

l> 

’.O 

( ) 

J 

30 

r» 

30 

0 


it*- PROGRAM TAILED TO RUM 


03 

2 . 7 1 

to 

3.04 

20 

4 . 53 

30 

9.90 

50 

K fc *f i? 

Ot > 3 

3. Oo 

ill 

'.09 

■35 

3 . 52 

3 0 

N -'£j 

2o 

DrXrk * 
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OR 


P(L> = N ( E ) (dE/dL) . 7.5 

THE KINETIC ENERGY E OF A PARTICLE OF MASS M IS GIVEN BY THE 
EXPRESSION, 

E = (GAMMA -1)*-Mi<-C 2 7.6 

WHERE 

GAMMA = SOR ( 1 / (1-BETA 2> ) ; BETA - V/C , 7.7 

AND C IS THE SPEED OF LIGHT. IT IS CONVENIENT TO EXPRESS ENERGY 
IN UNITS OF MC 2 SO THAT E - GAMMA-1. AFTER A LITTLE ALGEBRA WE 
FIND THAT 

E ~ L X/(X--1> T -• 1 , 7.0 

WHERE 

X -s (N 2)*( 1 - (L/T > > . 7.9 

US I NO FI TE SE DEE INI 1 1 ONS , 

dE/dL =- L (IM 2> / (2*1 ) I*! X*(X-t> 3 1 (-1/2) 7.10 

THE NORMALIZATION CONSTANT A IN EO 7.3 IS FOUND BY DIRECT 
INTEGRATION OF EO 7.3 TO BE 

A - (B-n/C LL ( B+l) ~ UL (-B+-1) 1 , 7. 3 1 

WHERE LL IS THE LOWER ENERGY LIMIT CONSIDERED FOR THE 
D3STR1BUT ION AMD UL IS THE UPPER LIMIT. PUTTING EOS 7.3, 7.8, 
7.10, AMD 7.11 TOGETHER, WE FIND AFTER A LITTLE ALGEBRA THAT 

P (L) " E (A>N 2 > / < 2 > 3 * C 1 / ( SOR ( X ) -SOR ( 1 -- X ) ) 3 R * 

30RU/X) > (X-t) (( B--3T/2) 7.12 

THIS DISTRIBUTION TN PULSE HEIGHT IS RATHER FLAT WHEN B=2.9, IS 
MONO TONIC UPWARD FOR B - 2.9 AND IS MONATOMIC DOWNWARD FOR B 
2.9. THE (B-3) EXPONENT TN THE LAST FACTOR OF EO >7.12 IS LARGELY 
RESPONSIBLE FOR THIS BEHAVIOIR . AN IMPORTANT FEATURE OF THESE 
CURVES IS THAT, AS B DECREASES FROM 3, THE STACI ING OF THE DATA 
POINTS JUST BELOW L (MAX ) BECOMES MORE PRONOUNCED. FOR B 3 , 

P(L) IS RATHER EVENLY MONATOMIC DOWNWARD FROM BEGINNING TO END. 

THE PROBABILITY OF ONE COUNT WITH L(I> < L ' L(H-1> IS 
P ( L ( I ) ) SO THAT THE PROBABILITY OF N(I> COUNTS IN THE INTERVAL IS 
P ( L ( I ) ) N ( I ) . A GOOD MAXIMUM LIfELIIIOOD FUNCTION P IS GIVEN AS 

P -= PRODUCT ( I ) P ( L ( I ) ) N ( I ) , 7.13 
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! <)h r 


ohml t > N(j ) i iju i r a ( i :• > j 


1-1 


nus rui*!(j] lon in easily max lmized by the: methods levh nrF.L' rn the 

AT s ACHED CUMPU ! hfR PROGRAMS. IF WE DO HOT BIN 1 1 10 DATA, A 
001 iVL N1 r:iMT RAX I HUH L j i LL I HOOD FUNOTiOH 15 GIVEN AS 


()L 


:;uh < ! ) lug pi i. >. j > i 


/„ t;. 


IT WE LET D AMD L ( MAX ) DC VARIABLES, WHERE L (MAX) MUC T BE 

ON AiT! v 1 1 IAN AIJY I 'H! SL HF.tOHf LU) USED, WE CAM LIMB B AND UMAX) 

wiu;.,i i maxim rzr cog r rn e the data, any answer 2 f< 4 is 

E'OSS : D'i E FOR SMALL VARIATIONS IN L (MAX ) > THE BEST AMD MOS'I 
BE L I LV7-.BL C RrSULTH WERT 00 I AIMED BY MAXIMIZING II LOG ND J /N WITH B 
AND I.HIAX) AO VmRIdDIFS, WHERE M IS l HE TOTAL MHMbER OT PARFLCIEfO 
U.T.D IN I III' MAXIMIZATION. Ti OHOUI I) BE NO IT D I MAT PiL) IS Rr, I III- R 
I I f . I AMD HIM AN’. I I INC I J I Tl WIN 01 1 1 1 AMD! )HLY El lANl II '0 Mil 0 f /I- 01 1 1 II. 

! JO! If f’UI 01.0 WILL NOT UIANOE ItJ ANY AFT RFC 3 ABl T IX MINI 1 1 IF SHAPE 

Ml E <1. ) PHI DM: II D BY EUHA1T0N 7.1*0 LXCEl'l A3 1 1 !L ENDS. DATA Wtl'll 
HD L (MAX ! IS RLI i.ECrr.D 1 1 IROUOi I L.(MAX) TO A LOWER ENERGY 
.'•EH :AX ) LIT). 1 1 II t ’.ROt )M! M T IN T AVOI i 01 1 1 1 i 0 if, AS f OU .UWS . I E WE 

r M. 1 ! 0._ ItJ Mil RJOlir EE P(l.) A MIRROR IMA OF P ' (1 > , I III. 

PRO .ViIML IMS’ I XCIliiNOi.D Bl I WEEN NIL TWO AT S I RJ Bi J I TONG ARE 
I II . HI SI- MIRROR' I liAfiES. SINCE EX Li I ANGLO ARE P.OUAL , TUT I.TAE AOT 0 

Al III.. Rl JlJi IDAhl 110 ALL OWN I01IMI 30 P'RLOLRVF I ' (I. ) II ! I iL" 

AR'i iil I I I C ! H) I'AEI III. 

ONE R'l S’O! 10 01 FEE .'HI MO 1 1 ARSON AMD IRON COSMIC KAYO BY 
'll IE ABOv'E MAX [MUM LtEM.TIIOOD 1RLATMLMT ARE SHOWN IN FJO 7. J 
■MID 'All. I 7.1. 

0 OWMMARY 


01 VI'EAL I )OEl UE II 01 IN I 001 0 WERE DEVELOPED 10 f X TRACT 
I 'ARAML I LEO I ROM DATA Oil RANDOM IVL-IMS. MIL BLOT UE 1 1 1!' II Wl KL 
BAOLD ON MAaIMOM l HI. I IIIOOD 0 ! A I LHE.I ITO. 1 1 IEOL ME. I HODS BEHOVED 
ONCOL SOI Ml JM I ELATING DA I A FROM LOOM JO RAYS FOR MIL' CASTS OF 
CAR BUM AND JE-'UN NUCIE.l. 
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TABLE 7.3 


ANALYSIS BF-DAT'A/iOtL'TRON NUCLEI 

A « J * • - * r ^ 


NUMBER OF DATA POINTS 82 

START I NO SPECTRAL INDEX FOR 1 HE MAX f MUM LH EL I HOOD ITERATION 
STARTING L (MAX ) 2114 

FINAL SPECIRAL INDEX 2. 09 ’8 
F [MAI L (MAX) 2022.0 
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I 


10 
20 
30 
40 
50 
60 
70 
80 
70 
1 00 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
750 
260 
270 
780 
, 70 
3 00 
3 j 0 
.'20 
3 30 
340 
350 

3 70 
"80 
390 
400 

4 10 
420 
4 30 
44 0 
450 
460 
4 70 
480 
49o 
500 
510 
520 
530 
540 
550 
560 
580 
590 
600 
610 
620 
630 
640 
650 
660 
6 •'0 


P ( L ) TO IRON DATA 
R. D. SHELTON 29 


' PROGRAM TO FIT 
' B: NASA48D. BAS 
CL3 

DIM L < 200 > , LL ( 200 ) , P# < 200 > , E# ( 200 ) 


UL=L (MAX ) - 1 L ( I ) 
JULY 85 09:32 


L (MAX ) REFLECTED 


DATA 1230,767, 1937, 1950, 1986, 1628,2466, 1367,2079,2243 
DATA 1974, 1895,2107, 1963,2195,1965,460,2023, 1602, 1958 
DATA 632, 1217, 1013,2034,1494,911, 1207,1948,1880, 1947 
DATA 1828,2003,944,901 , 1990, 1965, 1805, 1070,808, 1463 
DATA 1647,1370,625,640,1419,484, 1951,662,613 
DATA 3663, 1264,1330,2103,1726,1 125,2324,630,2468,1 133 
DATA 642, 1822,2176, 1953, 1154.754, 1841 , 1738,2259,977 


2034, 1874,2438, 1448,812,2241 


MOD 


- 0 THEN INPUT ZT 


DATA 3627,1533,1226,1730,1 
DATA 1590,952,2063 
3UM#=0 

F' 4 r-= "#####. : PP *•="### 

FOR 1=1 TO 82 

READ LL ( I > 

' SUM4=3UM# H_L# ( I > 

' PRINT I , LL ( I > : IF I 
NEXT I 

FOR 1=1 TO 82 

FOR J = I TO 82 

IF LL ( I ) LL < J ) THEM SWAP LL ( I > , LL ( J > 

NEXT J 

IF 3 MOD 3n =.». 

NEXT I 

' FOR 1=1 TO 02 

LF'RINT 1?" ’ 

' NL x r 1 

• PREFARL FOR 
CL.5 : PR I NT 
PRINT 


ORIGINAL PAGE IS 
OF POOR QUALITY 


■ THEN PRINT I: 

1 ; LL < I ) 

HERAT ION 


I1M— 3 s 
PRINT 

1 NPU I 

j 

PRINT 
IF Zt 

pr in r 

IF Zf 
fh IT 
‘ GQSUD 
3 l 

' 1 TERAT I ON 
BDSUB 810 
GOSUB 940 
LOO! -IT 
IF LOO! 

F'Rl NT- 
INPUT 


f=l: IT =- 1 r l-'RINT : INPUT 


ENTER STARTING EXPONENT 


, MU 


" ENT Eh. L (MAX " , Df! 
001 15: ET2#=-ET# 2 
INPUT " ENTER 3TE r ‘ FOR 
"" THEN HM4--VAL ( Z t ) 
INPUT " ETHER STEP FOF; 
. THEM HB4---VAL (2T) 
i THEN GOTO 460 


EXPONENT 


(MAX) 


7 =■ 


,zt 


IS : PRINT 


A 1 


XV 


SCHEME 


MOD 1 O : r 
U THEN 


'R I NT 
GOTO 


460 ELSE GOTO 500 


CR-GO ; C--CHANGC STEP 


"CHOICE 
PRINT 

IF Z-T“”00" THEN PRINT 1;’ 
IF ZT="C" THEN GOTO 380 
IF THEN GOTO 970 

IF Zr="00" THEN GOTO 600 
rr ZT=-"00" THEN INPUT Z-f 
GO ro 460 


IT;" EXP 


P-LPR1NT LL ( I > , L ( I ) ; 00-0UIT 

L (MAX ) = 
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z$ 


;B# 


END 
' Sfc:R 


TO COMPUTE SUMS FOR ITERATION 


G-S==0: LL=4f<S: 3 # 
~0R JJ = 1 TO 82 
[F LL ( JJ > - 
NEXT JJ 
XF-=E72*M ( j 
X 1 HET 24 H ( 1 


=Bft/ (1- 1 / ET 2# ) 


=B3*-.5 THEN L(JJ>=LL(JJ) ELSE L ( J J ) =2HB#-LL ( J J ) 


f-LL/L#) : ELL*=SOR(X#/ (X#-l) ) -1 
CD#- .5) / K # ) :EUL»=SOR(X#/ (X#-l) 


«' ; 1 /' * • * i i 1 J « L 1 ' 

f r 1 k " . I-L- 


Li- 1 . i it * \ j ; t r_ it ■ i :* 5 , ,» , 

*_/i — 4 . u j :! I Slit L; i i 4 

U ■ “ i l « •} r i i! , t« I - r r*« " ' 

f 1 i{ t ■ Ia ji j li <■ i. . II 1 t it— — l it i] t \ * I— I f u l t i 1 i A i 


) s » 
u w' ; 
-a - a - ul ‘ ( a 


1 1 i 

!• :! : 


i M i— ! 4 r> "f r _ lL 1 • i • •" K 1 


/ ! ,bt LU-L 


ir- 


"iT - , , *-f *• ‘ L- 


l/’> 

. ill 


I / i 

n J , } *■ i xil ' 

. ♦ “ • i i 

I » I { 


, - j ,, l * 1 "i * „ i\ i! , l. ,, „ ; 


i"; n • i i L? ’ i u 

• >t 1 Ml is". M.“ 


lLj kl i }-' i : 'J , 1 .‘Lt 1 


. tJ J I IU 


I 1 l* » * u 1 

w/U* fu 

i :: 1 r r 4 • 

of '/■>«;/* % 1“ 


{ i it / < 


i I l_i /-» i L. bi; i ■ { i a—' ; 1 .U ('-.lull M . • ; o 

-hiii i r-fkik-K,:.. . :i 

• I* 5 i:_ i s J f'ik = } ii'i 

71 IL.. * i .Mil !ii-fp '? 

*? ii i In'! - i-'i, hM 1 ' Hrih 'Zi Ji L'i! 


jj ; 


; » ; . v .{. 

b 1 »t“t {'■J ' z * ~ 

/f i>) «*; i . 

*. , ' i •)’ - 


'* '{ 


j rf .1 


i ; III !\ P -• [' J 
'ii . s; 

> *_{ (.*• s :-:i a - >. i-t 


1 n 

u! ! ' !i u -I 

:'! ! ! ! i M .JU . J i'fi |-;,T i >< .! 1 

L '"Z. i'. 'Pi'll UMKfi : i'-i; i I :: s i 'K I f : j l/.'l-.’n 
ivj r-T't ; i'1i I: H.' 


5 Pi; I!, i ,• :!V ! AT 


* 'DW 

!- 1 ! 

f )! \i } 

, •>-? , 

S OU 

1 

J “ 

c>{ , L . 

ill; 

let i 

97 ''‘ 

U.I/t 

i ’ 

1 f»t Jt- 

1 t i! 1 

, 0 


- ] 


i, : > 


ORIGINAL PAGE IS 
OF POOR QUALITY 
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